Tailoring the magnetic properties of Fe asymmetric nanodots 
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Asymmetric dots as a function of their geometry have been investigated using three-dimensional 
(3D) object oriented micromagnetic framework (OOMMF) code. The effect of shape asymmetry 
of the disk on coercivity and remanence is studied. Angular dependence of the remanence and 
coercivity is also addressed. Asymmetric dots are found to reverse their magnetization by nucle- 
ation and propagation of a vortex, when the field is applied parallel to the direction of asymmetry. 
However, complex reversal modes appear when the angle at which the external field is applied is 
varied, leading to a non monotonic behavior of the coercivity and remanence. 



I. INTRODUCTION 

During the last decade a great deal of attention has 
been focused on the study of regular arrays of magnetic 
particles produced by nano-imprint lithography. Besides 
the basic scientific interest in the magnetic properties of 
these systems, lithographed magnetic nanostructures are 
good candidates for the production of new magnetic de- 
vices, magnetic sensors and logical devices^ as high den- 
sity storage medial The properties exhibited by these 
nanostructures are strongly dependent on the geome- 
try, and therefore great control of the shape is funda- 
mental for the understanding and applications of such 
materials^. 

In the case of nanodots, two main reversal mechanisms 
for the magnetization have been observed, vortex nucle- 
ation and coherent rotatiorP. Vortex states are charac- 
terized by in-plane and out-of-plane magnetization. In- 
plane magnetization forms the vortex chirality, a clock- 
wise or counterclockwise magnetization rotation around 
a core. Out-of-plane magnetization defines the polar- 
ity, the up or down direction of the vortex core mag- 
netization. In this way, vortices exhibit four different 
magnetic states and can then store the information of 
four magnetic bits. Switching of the vortex polarity by 
the application of small magnetic fields along the dot 
axis suggests the possibility of using the vortex for high 
density storage media. Methods to control chirality in 
single FM layer elements exploit an asymmetry in the 
applied field, such as using a magnetic force microscope 
tipSlD ma g ne tic pulses^, or a magnetic field gradiemP, 
as well as the magnetization history^. Also the geome- 
try can be used to control polarity an d/or chirality, and 
dots with slight geometric asy mmet r y! 1 1 1 1 1 or triangular 
nanodot d 6 * 12 ! have been used for this purpose. 

On the other hand, some groups adopted the idea of 
asymmetric disks to achieve control over the vortex chi- 
rality with an in-plane magnetic fielcP^HIIl. Following this 
idea, Wu et aZPl studied vortex nucleation, annihilation 
and field distribution switching in 40-nm-thick NisoFe2o 
disk arrays, with a diameter of 300 nm and different de- 
grees of asymmetry. Their measurements and micromag- 



netic simulations showed that the nucleation and anni- 
hilation of vortices have a linear relation with respect to 
the ratio of the long/short asymmetry axes, while the 
switching field distribution oscillates. More recently, Du- 
mas et reported the synthesis and magnetic char- 
acterization of polycrystalline arrays of asymmetric Co 
dots where the circular shape of all the dots had been 
broken in the same way. In these arrays they showed 
how the vortices can be manipulated to annihilate at par- 
ticular sites under certain field orientations and cycling 
sequences. 

On another side, several works on magnetic nanopar- 
ticles focus on coercivity, which is strongly dependent on 
geometric parameters. Although current experimental 
facilities allow the fabrication of dots with a variety of 
geometries, precision is still a problem. Therefore, the 
possibility of controlling the coercivity by other means 
is highly desirable. Since coercivity is directly related to 
the reversal mechanism, one alternative is to induce dif- 
ferent reversal modes by modifying external parameters, 
such as the direction of the applied field. 

Following these ideas, in this paper, micromagnetic 
simulations have been used to investigate the angular 
dependence of the hysteresis and reversal modes for non- 
interacting asymmetric dots as a function of their ge- 
ometry. We focus on the behavior of the coercive field 
and remanence and conclude that magnetic fields applied 
along different directions on asymmetric dots are a pos- 
sible way of tailoring magnetic properties of nanodots. 



II. MICROMAGNETIC BACKGROUND 

The theory of micromagnetism was developed by 
Brown Jr. about 50 years agcE'. According to this 
model, a ferromagnetic system consisting of a large num- 
ber of individual magnetic spins is described using con- 
tinuous functions for the magnetization, the fields and 
the energies. Moreover, the amplitude of the magnetiza- 
tion vector M(r) has to be constant, but its orientation 
may change from one position to another. In this ap- 
proach, for a given magnetization distribution M(r) the 
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Gibbs free energy is 



G (M) = J dV 

2 1 

A (Vm) - ~HoM a (m • H dem ) - /Lt M s (m • H) 
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where H is an applied field, A is the exchange constant 
and Hdem is the demagnetizing field. Using the varia- 
tional principle to minimize the Gibbs free energy with 
respect to the magnetization M, the equation of the sta- 
ble equilibrium state is m x H e // = 0, where the effec- 
tive field H e // is defined as H e // = —dG/dM.. Because 
the magnetostatic interaction is of long range, analyti- 
cal solutions of micromagnetic problems can only be ob- 
tained for samples of simple shape and making use of 
simplifying assumptions^. Because of their symmetry, 
ring geome tries a re particularly suited for such analytical 
calculations^^. However, for systems having a complex 
geometry, like asymmetric rings, numerical simulations 
are required. 

The magnetization dynamics is governed by the 
Landau-Lifshitz-Gilbert equation (LLG)P^ 



dM 
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where 7 is the gyromagnetic ratio of the free electron 
spin and a is a phenomenological damping constant. 
The equation describes a combined precession and re- 
laxation motion of the magnetization in an effective field 
H e //. The calculations presented here have been com- 
puted by Oriented Micro Magnetic Framework software 
(OOMMFpB. Thus, the ferromagnetic system is spa- 
tially divided into cubic cells and within each cell the 
magnetization is assumed to be uniform. In order to as- 
sure a good description of the magnetization details, the 
size of the mesh has to be smaller than the exchange 
length of the material, defined by l ex = \/2A/ '/j-qM^ . 
The choice of the discretization scheme is validated by 
the fact that the numerical roughness (generated by the 
square mesh representation) corresponds to the real im- 
perfections on the lateral ring surface, arising for example 
from the resolution of the patterning methods used. 



III. SAMPLE SPECIFICATION 

Our starting point is a uniform circular dot with di- 
ameter d and height h. We introduce asymmetries in 
these dots by cutting specific sections characterized by 
the parameter S — r/R, as illustrated in Fig. 1. A sym- 
metric dot is characterized by 5 — 1 .0, while a semicir- 
cular one is represented by 5 = 0.0. To simulate the 
magnetic properties we used micromagnetic simulations, 
assuming that the dot-dot distance is large enough to 
consider every dot as independent^. For our simulations 



we use the typical Fe parameters: saturation magnetiza- 
tion M s = 1.859 x 10 6 A/m, exchange stiffness constant 
A = 45.78 x 10 -12 J/m, and a mesh size of 2 nm, where 
spins are free to rotate in three dimensions. Since we are 
interested in polycrystalline samples, anisotropy is very 
small and can be safely neglected^. In all the cases the 
damping parameter was chosen as 0.5. 




FIG. 1. (Color online) Geometrical parameters of a nanodot. 
The white surface represents the section that has been cut. 



IV. RESULTS AND DISCUSSION 

A. Magnetic field applied parallel to the x axis 

Our main concern in this work is to investigate the role 
of the asymmetric shape of the disk in the coercivity and 
reversal mechanisms of the magnetization. The hystere- 
sis curves for different diameters and 5, with h = 20 nm, 
are illustrated in Fig. 2, where we see for all diameters 
an almost square loop with coercivities strongly depen- 
dent on the geometry. For all diameters and 6 = 1.0 
the dots exhibit a small coercivity, which increases while 
decreasing S. 

In Figs. 3 and 4 we illustrate the general behavior 
of a dot for different values of h and 8. In these fig- 
ures we observe a non-monotonic behavior of coercivity 
and remanence which is the result of a competition be- 
tween exchange, local dipolar interactions and geometry 
in the region where the cut is made. Once we have a 
small asymmetry (<5 = 0.9) in some cases a decrease of 
the coercivity and remanence is seen. However, for fur- 
ther decreases of 5, remanence and coercivity are almost 
constant or grow continuously. 

In order to understand this behavior we look at snap- 
shots of the hysteresis and observe that the magnetiza- 
tion reversal occurs by a C formation followed by vortex 
nucleation and propagation. This behavior is illustrated 
in Fig. 5(a) for d = 60 nm, ft, = 20 nm and each of the 
three values of 5 used in Fig. 2. Looking at these snap- 
shots we can conclude that a vortex nucleates during the 
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FIG. 3. (Color online) Coercivity of asymmetric dots for dif- 
ferent height h, diameter d and S values. 



FIG. 2. (Color online) Hysteresis loops for asymmetric dots 
with height h — 20 nm for different diameters and 5 values. 



reversal for any value of 5. In symmetric dots, square 
loops are a sign of coherent reversal, and the appearance 
of a neck indicates that the reversal is driven by a vortex 
nucleation and propagation^. However for asymmetric 
dots reversal by vortex nucleation may lead to a square 
loop. 

Once a cut is made on the dot, a magnetic pole is 
formed on the new surface. Due to the Pole Avoidance 
Principle^, a C state nucleates in order to avoid the for- 
mation of a magnetic pole, followed by the nucleation of a 
vortex. However, when the asymmetry increases by cut- 
ting larger sections of the dot, the local effects described 
above compete with the new geometry, that is, the lack 
of circular symmetry, making vortex formation more dif- 
ficult. Then, in this last case, C formation and vortex 
nucleation become more difficult compared with the sit- 
uation when a small cut is made or a symmetric dot is 
considered, and an increase in coercivity and remanence 
appears. 



B. Angular dependence of the coercivity 

We also investigated the angular dependence of the 
magnetization, applying a magnetic field along a direc- 
tion defined by 9, which is the angle between the applied 
field direction and the x axis, as illustrated in Figure 1. 
Our results for d = 60 nm and ft, = 20 nm are depicted 
in Figs. 6 and 7. While symmetric dots show no angu- 
lar dependence, as expected, asymmetric dots exhibit a 
strong angular dependence. By increasing 9 the coerciv- 
ity decreases drastically until zero and a non hysteretic 
behavior is observed for 9 — 90°. This behavior allows us 
to conclude that shape anisotropy of an asymmetric dot 
may induce a hard axis of magnetization when 9 = 90°. 

Besides, the magnetization inside an asymmetric dot 
is oriented parallel to the x axis (easy axis), due to the 
shape anisotropy. When the applied field is reduced to 
zero, at remanence, each dot presents its magnetization 
along the axis, but it is measured at an angle 9 with 
respect to the easy axis. Then, one can approximate the 
remanence of an asymmetric dot by Mr(9)=Mrcos 2 (9), 
with Mr — Mr(9 = 0) the remanence measured at 9 
= 0. This behavior is an indication that asymmetric 
dots follow the behavior of uniaxial systems. Finally it 
is important to note that results for the remanence for S 
= 0.0 and 0.5 almost coincide. 
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FIG. 4. (Color online) Remanence of asymmetric magnetic 
dots as a function of 5. 

To understand this behavior, snapshots of the reversal 
of the magnetization are depicted in Fig. 5(b) and 5(c). 
It is seen that for S = 0.5 the reversal occurs by means of 
the nucleation and propagation of a vortex. For 0° < 8 < 
30°, the vortex nucleates at the center of the horizontal 
surface, and propagates along the y direction until it is 
annihilated at the opposite end of the dot. For 30° < 9 < 
60°, the vortex nucleates at the center of the horizontal 
surface and propagates along the x direction until it is 
annihilated at one end of the horizontal surface. For 
60° < 9 < 90°, a coherent reversal occurs. 

For S — 0.0 the reversal of the magnetization occurs 
through the nucleation of a C state and further nucle- 
ation and propagation of a vortex. For 0° < 9 < 60°, two 
regions of different magnetic orientation appear. These 
regions start growing until the full inversion of the mag- 
netization occurs. For 60° < 9 < 90°, a coherent reversal 
occurs. The non monotonic behavior of the coercivity 
in Fig. 7 can be explained by the existence of different 
reversal modes as a function of 9. 



V. CONCLUSIONS 

The results presented above show that the existence of 




FIG. 5. Snapshots of the magnetization of the dot during the 
reversal for a) 5=1.0, S — 0.5 and S — 0.0 when 9= 0, b) 
9=30, 60 and 90 when S = 0.5, and c) 6>=30, 60 and 90 when 
8 = 0.0. The arrows illustrate the direction of the magnetic 
moments. 



dot. In this case the coercivity and remanence are dras- 
tically modified as a function of 5 when a magnetic field 
is applied parallel to the x axis. Also we have extended 
our results to the case of an angular depedence of the co- 
ercivity and remanence, where a transition from vortex- 
mode to coherent-mode rotation has been observed. In 
this way asymmetry can be useful for tailoring specific 
magnetic characteristics of these systems. However, ex- 
perimental work remains to be done in order to observe 
this transition. 
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with diameter d — 60 nm and height ft = 20 nm for different 
angles and S values. 
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FIG. 7. (Color online) Coercivity and remanence of asym- 
metric dots with diameter d = 60 nm and height h = 20 nm 
for different angles and 5 values. 
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